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(N ■ ABSTRACT 

> ■ 

pf-v ' We have modelled the circumstellar envelope of o Ceti (Mira) using new observational con- 

V^ I straints. These are obtained from photospheric light scattered in near-IR vibrational-rotational 

CO ' lines of circumstellar CO molecules at 4.6 /im: absolute fluxes, the radial dependence of the 

Q^ I scattered intensity, and two line ratios. Further observational constraints are provided by ISO 

observations of far-IR emission lines from highly excited rotational states of the ground vibra- 
tional state of CO, and radio observations of lines from rotational levels of low excitation of CO. 
A code based on the Monte-Carlo technique is used to model the circumstellar line emission. 

The vibrational-rotational lines are sensitive to the radiation field, whereas the pure rota- 
tional lines, such as the rotational lines of low excitation measured at radio wavelengths and 
O I the rotational lines from highly excited states observed with ISO, are usually more sensitive to 

■4-J ' the temperature structure. These rotational lines have been the prime probe in most earlier 

j^ I investigations. 



We find that it is possible to model the radio and ISO fluxes, as well as the highly asymmetric 
radio-line profiles, reasonably well with a spherically symmetric and smooth stellar wind model, 
j^ ' However, it is not possible to reproduce the observed NIR line fluxes consistently with a 'standard 

H I model' of the stellar wind. This is probably due to incorrectly specified conditions of the inner 

regions of the wind model, since the stellar flux needs to be larger than what is obtained from 
the standard model at the point of scattering, i.e., the intermediate regions at approximately 
100— 400 R* (2"— 7") away from the star. Thus, the optical depth in the vibrational-rotational 
lines from the star to the point of scattering has to be decreased. This can be accomplished 
in several ways. For instance, the gas close to the star (within approximately 2") could be in 
such a form that light is able to pass through, either due to the medium being clumpy or by the 
matter being in radial structures (which, further out, developes into more smooth or shell-like 
structures). Further observations of the gas in the stellar wind close to Mira are required to 
resolve this problem. 

The model circumstellar envelope, which reproduces the observables reasonably well, has a 
mass-loss rate of 2.5x10~^Mq yr~^, and a turbulent velocity of l.Skms""'^, given a terminal 
expansion velocity of the wind of 2.5kms~^. 

Subject headings: stars: AGB and post-AGB — circumstellar matter — late-type — mass-loss — infrared: 
stars 



1. INTRODUCTION 

The inner regions of circumstcUar envelopes 
around asymptotic giant branch stars are still 
enigmatic. These regions, located further out 
than the pulsating outer atmosphere but typically 
within 100 stellar radii, are of interest for several 
reasons; The circumstellar wind is accelerated here 
to its terminal velocity, the chemistry is 'frozen 
out' approximately at its photospheric values, and 
dust nucleation and growth take place. This is 
important for the mass loss of the star, a pro- 
cess we do not understand how and why it occurs. 
For example, we do not know whether matter is 
ejected in puffs or in a steady, homogeneous wind. 
The mass loss returns matter to the interstellar 
medium and plays a vital role in the cosmic cy- 
cling and enrichment of matter, see for example 
Gustafsson and Ryde (2000). In order to study 
these inner regions we need to observe them in de- 
tail. However, very few observational constraints 
exist on the conditions in these complex regions. 

By combining observations, sampling gas lo- 
cated at a distance on the sky of approximately 
1" and beyond, and a detailed modelling of the 
circumstellar envelopes one might be able to put 
constraints on the inner regions. Ryde et al. 
(2000) observed scattered photospheric light in cir- 
cumstellar vibration- rotational lines of CO 2"— 7" 
away from Mira. They analysed the emission lines 
from the circumstellar envelope using a simple an- 
alytic approach, an analysis which suggested a 
time constant mass-loss rate, when averaged over 
100 years, over the past 1200 years. In this pa- 
per we present a detailed analysis both of these 
observations and of those which tell us about the 
properties of the whole envelope. The modelling 
presented here provides a better insight into the 
conditions close to the star, and gives a larger de- 
gree of flexibility to test different models. 

o Ceti (Mira or HD 14386) is the prototype of 
the mira class of long-period variables with char- 
acteristic pseudo-periodic brightness variations in 
time. They are cool, red giants on the asymptotic 
giant branch (AGB) with considerable mass-loss 
rates and circumstellar envelopes (CSEs) of large 
extension. The visual-light variations are related 
to large-amplitude radial pulsations which cause a 
large variation of the effective temperature. Note, 
that o Ceti is a binary system. The angular dis- 



tance between Mira A and the hot, compact com- 
panion star Mira B (VZ Ceti) is 0.6" (Karovska et 
al. 1993; 1997). The companion could be a white 
dwarf with a mass of about 1 Mq , a luminosity 
of 2 Lq, and a temperature of more than 30000 K 
[these parameters are based on an analysis of In- 
ternational Ultraviolet Explorer (lUE) spectra by 
Stickland et al. (1982)]. Danchi et al. (1994) sug- 
gest it is embedded in an accretion disk, giving 
rise to an abnormal illumination of Mira A. 

In general, the kinematics, the density and the 
temperature structures of CSEs of AGB stars seem 
to be relatively simple, with an over-all spheri- 
cal symmetry prevailing (see, e.g., Olofsson 1996). 
The stable and abundant CO molecule is self- 
shielding (Glassgold 1996) and therefore of ap- 
proximately constant abundance trough out the 
envelope. However, the chemical structure is cer- 
tainly quite complex, restricting us to schematic 
modelling of these objects [see for example Glass- 
gold (1996)]. The lack of spatial resolution in the 
observations, nonetheless suggesting more com- 
plex geometries (e.g., Stanek et al. 1995), also 
hampers the analysis. 

The observations of CO lines at millimeter 
wavelength in Mira's wind indicate some com- 
plexity. An asymmetry is clearly visible in the 
radio- line profiles. A number of suggestions in 
the literature have been made for multi-wind sce- 
narios, such as winds of different characteristics, 
jets, or other phenomena which introduce dra- 
matic variations in the envelope structure with 
distance. Thus, several combinations exist of 
mass-loss rates, expansion and turbulent veloci- 
ties, which are able to reproduce the radio-line 
profiles fairly well. This means, for instance, 
that there is no consensus as regards the ac- 
tual expansion velocity of Mira's wind. Crosas 
et al. (1997) experiment with an expansion ve- 
locity (uc) of approximately 2kms^^ and a tur- 
bulent velocity {vt) of approximately 4kms^^ 
in their version of the inner wind, while Young 
(1995) arrived at Wo=4.8kms^^. However, Young 
(1995) points out that the wings show an expan- 
sion velocity of lOkms"^. Knapp et al. (1998) 
model a fast outer wind with a mass-loss rate of 
M=4.4x10~^Mq yr^^, and an expansion veloc- 
ity of 6.7±1.0kms~^. An inner wind is supposed 
to be a resumed wind in analogy with the de- 
tached shells found in four carbon stars (Olofsson 



et al. 1996). This slow wind component has a 
lower mass-loss rate and a much lower expansion 
velocity than what is found normally in Miras; 
M=9.4xlO~*^Moyr~i and Wc=2.4±0.4kms-i. 
Planesas et al. (1990a) interpret the CO J=l-^0 
and J =2—^1 lines as partly originating from a 
spherically symmetric CSE with Uc=3kms^^. 

Numerous high-resolution, optical and infrared 
wavelength measurements provide evidence that 
Mira itself is elongated (e.g., Karovska et al. 1991, 
1997; Haniff et al. 1992; Wilson et al. 1992; Quir- 
renbach et al. 1992). Karovska et al. (1997) sug- 
gest that the asymmetry could be due to unre- 
solved bright spots or non-radial pulsation, or even 
a wind interaction with Mira B. As an alternative 
to an asymmetric atmosphere, Danchi et al. (1994) 
discussed the possibility of the observed asymme- 
tries being caused by non-uniform dust shells close 
to the star. Also Planesas et al. (1990a, b) inter- 
preted their asymmetric CO J=1^0 and J=2-^l 
line profiles and brightness maps, which are at a 
spatial resolution of 6", as originating partly from 
a very slow, partially coUimated wind, due to high 
densities of gas in the equatorial plane close to the 
star. 

We will, nevertheless, assume a spherically sym- 
metric model for the CSE since the situation is 
unclear as regards the physical origin of the com- 
plexity, and its dependence on the distance from 
the star. Also, Ryde et al. (2000) concluded from 
their near-IR (NIR) CO line data that they do not 
necessarily indicate a non-symmetric, inner wind. 

2. OBSERVATIONAL DATA 

Our modelling of Mira's wind will be compared 
with three sets of observations probing different 
regions of the CSE. 

First, spatially and spectrally resolved, vibrational- 
rotational line emission from the CO fundamen- 
tal band lying at 4.6 /xm is used as a probe of 
the inner regions (Ryde et al. 2000). The lines 
observed are R(l), R(2), and R(3). These emis- 
sion lines consist of photospheric light scattered 
by CO molecules in the intermediate regions (ap- 
proximately 100-400 R,) of the CSE around Mira. 
These observations are the most important ones 
in the discussion presented here, providing a new 
constraint to the modelling of the CSE. The emis- 
sion lines were detected in observations with the 



Phoenix spectrometer mounted on the 4 m Mayall 
telescope at Kitt Peak. The reduction and ab- 
solute flux calibration of these observations are 
described by Ryde et al. (2000). The uncertain- 
ties in the absolute fluxes are probably lower than 
±50%, which we regard as a conservative estimate. 

Second, as a further constraint on our model we 
have retrieved far-IR (FIR) observations of Mira, 
which probe regions within 2" (< 100 R,) of the 
star, obtained with the Infrared Space Observa- 
tory (ISO) (Kessler et al. 1996). These high- lying 
rotational lines of ^^CO are retrieved from the 
ISO-data archive^. The Long Wavelength Spec- 
trometer (LWS, Clegg et al. 1996) was used in 
the grating mode (LWSOl) providing a mean spec- 
tral resolution element of approximately 0.7 /ini. 
We will discuss only the line-integrated fluxes. 
The 32 minutes of observations were performed in 
the morning the 4**^ July 1997 during ISO rev- 
olution 596. The reductions were made using 
the most recent pipeline, basic reduction package 
OLP (v. 7) and the ISO Spectral Analysis Pack- 
age (ISAP v.1.5). The pipeline processing of the 
data, such as the flux calibration, is described by 
Swinyard et al. (1996); the combined absolute and 
systematic uncertainties in the fluxes are large, 
of the order ±50%. Mira is considered to be a 
faint source in the ISO-LWS reduction procedures 
(Flux< 50 Jy at 120 /im). The observations con- 
sist of weak emission lines on a strong contin- 
uum. A large uncertainty is introduced in the 
subtraction of the continuum. We have measured 
the emission intensities in the same manner as is 
described in Ryde et al. (1999c) for IRAS 15194- 
5115. The lines studied are J=19^18, J=18^17, 
J==17^16, J=16^15, at 137.2, 144.7, 153.3, and 
162.9 /xm, respectively. The line-integrated fluxes 
measured yield a mean flux of approximately (2 ± 
l)xl0~^° Wcm~^ for the four hues. 

Third, radio observations, made with the James 
Clerk MaxweU Telescope (JCMT) on Hawaii, of 
the CO J=2->1, J=3-^2, and J=4^3 lines are 
retrieved from the public archive of the JCMT^. 
These data are taken at face value. By compar- 



^This project is partly based on observations with ISO, an 
ESA project with instruments funded by ESA Member 
States (especially the PI countries: France, Germany, the 
Netherlands and the United Kingdom) and with the par- 
ticipation of ISAS and NASA, http://isowww.cstcc.csa.nl/ 

^http://www .jach.hawaii.edu/ 



ing different observations we estimate the uncer- 
tainty in the intensity scale to be approximately 
±20%. The intensity scale is given in main-beam 
brightness temperature, Tmb- These data provide 
constraints on the outer regions of the CSE, i.e., 
approximately 1000 R, away from the star and be- 
yond. 

In Tables 2, and 3 the observational results of 
the NIR, ISO, and JCMT observations are pre- 
sented. In Table 2, the mean values of the NIR 
line intensities obtained west, east, north, and 
south (Ryde et al. 2000) are given, as well as the 
R(l)/R(3) and R(2)/R(3)-linc ratios. The inten- 
sities are means for a position (2±0.5)" away from 
the star, the ratios are means over the distance 
2—3.4". Table 2 also presents the radial bright- 
ness distribution of the NIR line emission given 
as the exponent of the intensity power-law (dlog 
I/dlog/3) of the R(2) vibrational-rotational line 
and is a mean over the entire range (/3 is the angu- 
lar distance from the star on the sky). The uncer- 
tainties quoted for the line ratios and the slope are 
pure measurement uncertainties. The R(l)/R(3) 
ratio in Table 2 may be uncertain due to system- 
atic uncertainties in R(l), cf. Ryde et al. (2000). 
The ISO flux in Table 3 is given as the mean of the 
four measured lines in order to reduce the uncer- 
tainty since the data is noisy and the line intensi- 
ties should not vary much. 

3. THE MODEL OF THE CIRCUM- 
STELLAR ENVELOPE 

We have used a Monte-Carlo (MC) approach to 
the multi-level, non-LTE, radiative-transfer prob- 
lem when modelling the circumstellar wind of 
Mira. We have applied the methodology of Bernes 
(1979) to the transfer of line-radiation of CO in a 
spherically symmetric and smooth CSE, formed 
by a constant mass-loss rate, and expanding at a 
constant velocity. The inner and outer radii of the 
CSE are used as the input for the modelling. For 
a more detailed description of the modelling, see 
Schoier (2000). 

The model is constructed by iteratively solving 
the equations of radiative transfer and the statis- 
tical equilibrium equations for the level popula- 
tions of ^^CO. The solution is found when a steady 
state of the populations, as a function of radial 
distance from the star, is achieved. The source 



function of the entire wind is thereby known, 
which makes possible a straight-forward calcula- 
tion of the spectrum by the solution of the radia- 
tive transfer through the envelope along the line- 
of-sight, resulting in predicted intensities of the 
observed spectral lines as a function of the dis- 
tance from the star. 

A draw-back of Monte-Carlo simulations in gen- 
eral is the slow convergence, typically scaling with 
\/Nitor- On the other hand, the advantages are 
several. The programming can be done in close 
relation to the physical reality; the radiation field 
can be taken into account in a correct manner, 
and inhomogeneities and complex geometries are 
relatively easily implemented. Velocity fields and 
gradients of any size are allowed in the program- 
ming. 

In our MC model, complete angular and fre- 
quency redistribution is assumed, the absorption 
and the emission profiles being the same. Thus, 
the local emission is assumed to be isotropic and 
the scattering is assumed to be incoherent. The 
distribution of frequency shifts given by the emis- 
sion profile is defined by a gaussian Doppler pro- 
file: 



*M = 



: CXp 



— (^ — ^'o — u • nvo/c)"^ 



(1) 



The width of this profile, ct, is given by the lo- 
cal, thermal broadening and the micro-turbulent 
broadening. The velocity field in the expanding 
wind will allow shifts from the line centre as given 
in the laboratory frame, fp. The velocity field is 
represented by the vector v and the direction of 
the beam, i.e. the random direction of one emit- 
ted photon, by the unit vector n. 

For the rotational transitions, with spontaneous 
de-excitation rates of the order of 10~'^s~"'^, the 
assumption of complete frequency redistribution 
(CRD) is a good assumption, since the collisions 
will have time to reset the frequency-memory of 
the excited atoms. However, for the vibrational- 
rotational transitions, with rates of the order of 
10 s^^, this assumption may be a fallacy. For these 
infrared transitions, photons are more likely coher- 
ently scattered in frequency. The more scattering 
events, i.e., the higher the radial optical depth in 
the wind, the more important it is that the model 
is made to properly take into account coherent 



scattering in the model. 

By dividing every NIR transition into a large 
number of subdivisions, but conserving the over- 
all line profile, we are able to simulate coherent 
scattering in the MC code. The absorption and 
the emission profiles are assumed to be the same. 
We find that for optical depths of the order of 
50 or less, the difference between a CRD and a 
coherent description of the NIR transitions is less 
than the numerical noise, which is intrinsic to the 
MC method. For a Mira model with a CSE with 
an inner radius at 13 R, , the radial optical depth is 
approximately 40, which motivates our approach 
to run the model with all transitions in CRD. 

4. MODEL PARAMETERS 

Measurements of the multi- wavelength angular 
diameter of Mira (Haniff et al. 1995), combined 
with trigonometric parallaxes measured with the 
Hipparcos satellite (Van Leeuwen et al. 1997), sug- 
gest a mean radius of (464±80) R© . The distance 
to Mira is (UStf^) pc (ESA 1997). iFrom Mahler 
et al. (1997) one can derive a bolometric lumi- 
nosity Ltot=8900LQ, a temperature Toff=2400K 
and a radius of 55OR0=3.8xlO^'^cm for the phase 
of o Ceti at October 29, 1998, the date the NIR 
observations were made, see Figure 1. Mira was 
at its maximum radius at this epoque according 
to Mahler et al. (1997). The ISO and radio data 
were observed at other times, but are not as sen- 
sitive to the radiation field, and thereby to the 
phases of Mira, as the NIR data. These val- 
ues of the model parameters of Mira are consis- 
tent with the corresponding values obtained by, 
among others, Danchi et al. (1994) and Haniff 
et al. (1995). These authors arrive at similar tem- 
peratures; a mean (over phases) effective tempera- 
ture of rcff^2800 K. The terminal expansion veloc- 
ity is set to be Wo=2.5kms~^ (Knapp et al. 1998). 
The model parameters we have adopted are given 
in Table 1. 

The radiative energy input to the wind system 
originates from the central star and is specified 
by a luminosity and an effective temperature of 
the star and by a luminosity and effective tem- 
perature of the surrounding dust, cf. Tablet. 
Since the dust is concentrated close to the star 
(Bester et al. 1991), the ambient radiation from 
the dust further out in the CSE is considered to 



be small and it is not taken into account in our 
model. The dust will absorb part of the stellar flux 
and re-radiate it in the infrared as thermal radia- 
tion. Lopez et al. (1997) fit their observations with 
two spherically symmetric dust shells at 3 and 10 
stellar radii, with inhomogeneities or clumps em- 
bedded in the spherically symmetric dust shells. 
Danchi et al. (1994) fit their observations with a 
dust distribution which peaks at 2.7 stellar radii 
away from the star, and where the temperature 
is 1050— 1300 K, depending on phase. Also Bester 
et al. (1991) deduced an inner extension of the dust 
shell of about three stellar radii, and a dust tem- 
perature of 1200 K. Thus, the spatial dust distri- 
bution peaks at approximately 3 stellar radii with 
a temperature of about HOOK, leading to an up- 
per limit of the re-radiated luminosity of 30% of 
the total bolometric luminosity; 2700 L0. Haniff 
et al. (1995) found, however, that the circumstcl- 
lar dust in Mira contributes less than 20% of the 
bolometric luminosity. We therefore set the dust 
luminosity to be 1800 Lq, implying a stellar con- 
tribution of 7100 Lq. The former radiation, which 
is mainly in the infrared, is important for pump- 
ing vibrational states, leading to non-thermal level 
populations of the rotational states which could 
affect the NIR line intensities. Hence, we model 
the inner energy source with a spectral energy dis- 
tribution consisting of two black-bodies, using the 
temperatures and luminosities above, see also Ta- 
ble 1. 

The kinetic temperature structure of the gas af- 
fects the CO level populations through collisions, 
modelled through coUisional rates [see the discus- 
sion by Ryde et al. (1999c)] which are functions of 
the temperature. The temperature also affects the 
thermal broadening of the absorption and emission 
profiles. These are, however, mostly determined 
by the larger turbulent velocity. The gas temper- 
ature structure of the CSE can be calculated in 
a self-consistent way in our model, by solving the 
energy balance equation of the system (cf. Ryde et 
al. 1999c). However, in the present discussion we 
select, for simplicity, a structure of the gas kinetic 
temperature in the wind according to: 



T{r) = n 



/lO^'^cm 



C 



(2) 



where the exponent C~li see for example Kahane 
and Jura (1994). Tq is the temperature at 10^^ cm 



and r is the distance from the star in centimeters. 
This kinetic temperature structure and the exci- 
tation temperature structures calculated by the 
model are shown in Figure 3b. 

Mira B certainly also adds mechanical energy 
to the circumstellar environment, which will lead 
to an additional heating. The NIR lines are, how- 
ever, radiatively excited and not greatly affected 
by the kinetic temperature of the surrounding 
gas/dust. The radio lines may, however, in prin- 
ciple be affected by their kinetic environment, but 
are formed further out where the influence of the 
binary star can be considered to be small. 

The CSE is characterised by the mass-loss rate 
of the star and the expansion velocity of the wind, 
which together give the density as a function of 
the radial distance from the star. The density of 
CO will depend on the fractional abundance of CO 
relative to the abundance of hydrogen molecules, 
H2. This fractional abundance of ^^CO molecules 
is assumed to be /co=5xlO~'' (Knapp et al. 1998) 
throughout the wind. This is a value generally 
used for oxygen- rich envelopes, and the same value 
was adopted in Ryde et al. (2000). In the case of 
Mira the mass-loss rate is relatively small, leading 
to an increased importance of the radiation field, 
especially for the NIR lines; see the discussion on 
the dusty envelope of the object IRAS 15194-5115 
with a high mass-loss rate, where the vibrational- 
rotational transitions arc found to be unimportant 
(Ryde et al. 1999c). 

The ^^CO molecule data are the same as those 
used by Ryde et al. (1999c). We use 30 rotational 
levels for each of the ground and first vibrational 
states. This is generally sufficient in order to treat 
the CO excitation in CSEs properly, cf. Ryde 
et al. (1999c). The colfisional rates of the NIR 
vibrational-rotational transitions are much smaller 
than the radiative de-excitation rates and they are 
set at zero in the model (see for example Ryde et 
al. 1999a). The colfisional rates of pure rotational 
transitions are non-zero. Kirby-Docken and Liu 
(1978) provide oscillator strengths, which for our 
transitions are consistent with Hurc and Roucff 
(1996)3. 

The photo-dissociation radius of the CO enve- 
lope (rp) is determined by the interstellar UV-field 
(Mamon et al. 1988). In the case of Mira, this 



^Note the missing cube in their Eq. 3 



value is relatively uncertain. We apply an outer 
radius of the CO shell of 2x 10^^ cm, which corre- 
sponds to a diameter of about 20" at the distance 
of Mira of 128 pc, cf. Planesas et al. (1990b). 

5. RESULTS PRODUCED BY THE MODEL 

For objects with low mass-loss rates, the NIR, 
vibrational-rotational lines are sensitive to the 
stellar intensity, whereas pure rotational lines, 
such as the low rotational lines measured at ra- 
dio wavelengths and the higher rotational lines ob- 
served with ISO are usually more sensitive to the 
temperature structure. The vibrational-rotational 
lines are therefore probes of the light-absorbing 
properties of the inner regions. 

The estimated fluxes given by the model are 
sensitive to many of the input parameters. We 
have run our model by fixing as many parameters 
as possible from the literature. Our variable pa- 
rameters are the temperature (Tq) at a distance 
from the star of 10^^ cm, the mass- loss rate (M), 
the terminal expansion velocity (we) and turbulent 
velocity (wt), as well as the inner radius of the CSE 

(n). 

Initially we apply a simple, 'standard model' of 
a spherically symmetric and smooth CSE expand- 
ing at a constant velocity. Our observational con- 
straints on the structure of the outer part of the 
wind are given mainly by the radio data, and to a 
lesser extent, the far-IR ISO data in the same man- 
ner as by Ryde et al. (1999c). We later confront 
this model with the new observational constraints 
consisting of the absolute fluxes, the radial depen- 
dence of the scattered intensity as a function of 
angular distance, and the two line ratios obtained 
from photospheric light scattered in vibrational- 
rotational lines of circumstellar CO. 

5.1. The 'standard' model 

Previous radiative transfer models analysing 
the CSEs around bright M stars (e.g., Kahane and 
Jura, 1994) assume large velocity gradients (LVG). 
The model presented here is a more general code, 
based on the MC method, with a detailed radiative 
transport, including non-local effects, and treating 
the gas in full non-LTE. In the limit where the ra- 
diation field does not have a large influence [see 
for example the discussion on the 'dusty' carbon 
star IRAS 15 194-5115 (Ryde et al. 1999c)], non- 



local effects are not very important. Kahane and 
Jura (1994) chose an inner radius of 5x10'^'* cm, 
i.e., approximately 13 R*, which is relatively far 
out in the wind. In our standard model we adopt 
this value. Note that the dust is supposed to be 
situated well within this distance. 

The turbulent velocity and the expansion ve- 
locity of the gas in the CSE are free parameters in 
the model which we tune together with the other 
free parameters in order to explain the observed 
asymmetric radio-line profiles by self-absorption 
in an optically thick wind, i.e., the mass- loss rate 
has to be sufficiently large. The temperature at 
10^'^ cm (see Eq. 2) we set to ro=150K [a value 
consistent with those derived from a large survey 
of optically bright carbon stars (Schoier & Olof- 
sson 2000) where the kinetic temperature struc- 
ture is obtained in a self-consistent manner] . The 
mass-loss rate is found to be 2.5x10^^Mq yr^^, 
the terminal expansion velocity is 2.5kms^^, and 
the turbulent velocity required is 1.5 km s""'^, see 
Table 1. 

In Figure 2 and 3 the model results are shown 
and compared with some of the observations. We 
find that it is possible with this model to meet 
the observational constraints provided by the ra- 
dio and ISO data reasonably well. The three ra- 
dio lines are reasonably well reproduced by the 
model (Figure 2a-c) , suggesting a certain impor- 
tance of self-absorption in the CSE around Mira. 
Also the modelled fluxes of the lines observed by 
ISO are reasonably well reproduced (Figure 2d). 
The scattered, NIR fluxes, as a function of the 
distance from the star, are shown in Figure 3a and 
they decline roughly as a power law, with a slope 
of —3.4 for the R(2) transition. Figure 3b shows 
the assumed kinetic-temperature law, as well as 
the excitation temperature for some selected lines, 
throughout the envelope. The excitation temper- 
atures of the radio lines show that the populations 
of the lower lying levels of the scattered infrared 
lines R(l), R(2), and R(3) are not in local ther- 
modynamic equilibrium (LTE). 

In Figure 3c the tangential optical depth is 
shown and gives an indication of the depths-of- 
formation of the different lines. Furthermore, in 
Fig. 3d the relative population of level i for the 
different rotational levels of the two lowest vibra- 
tional states of CO for our model arc shown as a 
function of the distance from the star. The popu- 



lations of the rotational states of the first excited 
vibrational state are several orders of magnitude 
lower than those of the ground level. Furthermore, 
it can be seen in the figure that these populations 
vary significantly with the distance from the star, 
a fact that can not be taken into account in the 
analytic approach of Ryde et al. (2000). 

The angular behaviour of the vibrational- 
rotational line-intensities and the line ratio of the 
two R-branch transitions are reasonably well re- 
produced (Table 2). However, the intensities come 
out significantly too weak, approximately a fac- 
tor of 5—10 lower than observed. The maximum 
optical depth in the line-of-sight is approximately 
25, while the radial optical depth is approximately 
a factor of two higher. We find that by varying 
the mass-loss rate, the temperature structure and 
the turbulent velocity, it is not possible to model 
consistently all the observations. 

By introducing a velocity law, i.e., an accel- 
erating wind, it is possible to increase the NIR 
fluxes derived from the model. Assuming a lin- 
early increasing, with radius, expansion velocity 
of the CSE, NIR fluxes that are about 2 times 
larger than in the standard model is possible to 
obtain. 

5.2. The 'cavity' model 

We find that it is impossible to reach the high 
fluxes measured in the NIR lines far out in the 
CSE with a standard wind model extending all 
the way to the vicinity of the star. There are nu- 
merous ways of explaining an observed NIR flux 
that is lower than expected, but for a higher flux 
there are not many alternatives. The key problem 
is to prevent the innermost region of a smooth en- 
velope to lose the radiation by scattering events 
in a geometrically diverging velocity field. In or- 
der to increase the NIR fluxes one can introduce 
an inner, empty region - i.e., postulate a large in- 
ner radius, which will increase the NIR flux at the 
point where the observations are made. 

By moving the inner radius of the CSE fur- 
ther out, high-density material close to the star 
is removed allowing the central radiation-fleld to 
penetrate further out into the wind. This will in- 
crease the excitation of the CO molecule, includ- 
ing the low-energy rotational transitions. With 
a cavity it is possible to explain reasonably well 



the FIR fluxes as well as the observations of the 
NIR vibrational-rotational lines of CO, cf. Ta- 
bles 2 and 3. However, the radio- line intensi- 
ties will come out to strong in the model, due 
to the efficient radiative excitation, and the self- 
absorption present in these lines will become less 
pronounced since the optical depth decreases (by 
approximately a factor of 5 in the radial direction). 
We find that an inner boundary at 2.5 x 10"'^^ cm 
(approximately 1" or 53 stellar radii) combined 
with a mass-loss rate of M=2.5xlO~^M0 yr~^, 
i.e., the same as in 'the standard model', will re- 
sult in a sufficiently high flux of stellar radiation 
reaching the wind regions beyond 2" and at the 
same time enough scatterers at the location of our 
observations, i.e., at angular distances of 2"— 7" 
away from the star. The ISO lines are formed fur- 
ther in than the NIR vibrational-rotational lines, 
and will start to become sensitive to the adopted 
inner radius. In the modelling, we have chosen 
the same temperature structure, velocity structure 
and turbulent velocity as was obtained from 'the 
standard model'. 

The results of our best 'cavity model' of the 
CSE around Mira are presented in Tables 2 and 3, 
where the modelled intensities, line intensity ra- 
tios, intensity slopes, the modelled radio-fluxes, 
and the modelled mean ISO-flux arc compared 
with the observations. 

The modelled intensities of the three NIR lines 
[R(l), R(2), and R(3)] all decline in a log/-log/3- 
plot with a slope of approximately —3, which is in 
good agreement with the observations. Note that 
there are systematic uncertainties in the R(l) line 
due to the detector, which will lead to an uncer- 
tain R(l)/R(3) ratio [see the discussion by Ryde 
et al. (2000); this has to be borne in mind when 
discussing the line ratios]. 

6. DISCUSSION 

The Kitt Peak observations provide our mod- 
elling with new, strong constraints. These arc the 
NIR absolute fluxes of the three lowest vibrational- 
rotational transitions in the R-branch of CO, their 
radial distributions, and the two line ratios from 
the three vibrational-rotational lines observed. 

We find that it is not possible to meet these new 
constraints with a 'standard model' of the wind, 
calculated on the basis of the several far-IR ro- 



tational transitions within the ground vibrational 
level of CO as observed by ISO, and radio obser- 
vations of millimeter rotational transitions of CO; 
our modelled NIR fluxes are an order of magnitude 
too low. 

By introducing a cavity in order to simulate a 
higher flux at the point of scattering and keeping 
the mass-loss rate relatively high, however, we are 
able to model the NIR fluxes and the brightness 
distribution, i.e., the decline in the intensity as a 
function of the distance from the star, providing 
a slope (dlog I/dlog/3) of approximately —3. This 
results in few scatterers on the way and hence in- 
crease the number of photons that will be able 
to be scattered, and, on the other hand, main- 
tains enough scatterers at the point of our obser- 
vations for the modelled fluxes to be large enough. 
The important free parameters that couple to each 
other to give the number density of scatterers are 
the inner radius (r;), the mass- loss rate (M), and 
the expansion velocity (we). 

Thus, in order to reproduce the observed, high 
NIR fluxes, the radial optical depth should not 
be too large. If there are too many scatterers on 
the way radially, many photons will scatter 'side- 
ways' due to the shifts out of the absorption pro- 
file caused by the velocity of the fiow at a point 
close to the star and will never make it out to 
our point of observation. The high optical depths 
trap the photons as was suggested already by Dyck 
et al. (1983). Thus, the velocity field is of great 
importance. If the geometrical velocity gradients 
are large across one optical depth unit (the mean 
distance between scattering events), the absorp- 
tion profiles of the absorbing molecules will be 
shifted with respect to the emission profile of emit- 
ting molecules, leading to a decrease in the optical 
length for the emitted photons. This will lead to 
a loss of photons and a steepening of the decline 
of the intensity with distance from the star. 

An important quantity for the fiux levels at a 
certain point in the wind is the relation between 
the expansion velocity, and the turbulent velocity 
(providing most of the width of the inherent line 
profiles). If the turbulent velocity, which is of the 
order of 1 kms~^, is comparable to the wind veloc- 
ity, the velocity fields will be very important. Also 
in the case of Mira, this makes the modelling of the 
CSE sensitive to the ratio between the turbulent 
velocity and the fiow velocity. This is also what 



we find; a sensitive parameter in the modelling is 
this ratio. 

A physical reason for a higher flux could be 
a void of CO molecules, for instance due to an 
excess UV flux from the white-dwarf companion- 
star that certainly will partly dissociate CO, or 
a cavity caused by a decline of the stellar wind 
to a lower expansion velocity and a lower mass- 
loss rate, resembling the double wind scenario of 
Knapp et al. (1998). Other reasons could be that 
the gas close to the star is in such a state that 
light is able to pass through relatively unaffected, 
either due to the medium being clumpy with a 
low filling factor, or by the matter being in ra- 
dial structures in the innermost part of the CSE, 
which, further-out, developes into more smooth or 
shell-like structures. These structures could occur 
depending on the expulsion mechanism of matter, 
such as a non-homogeneous mass-loss due to large 
convective eddies or magnetic spots, see for ex- 
ample Soker and Clayton (1999). Further reasons 
could be CO molecules absorbed onto dust grains, 
which are evaporated further out, or a small ra- 
dial velocity gradient (i.e., the velocity slowly ap- 
proaches the terminal velocity of the wind) . 

For large radial optical depths, the slope (dlog 
I/dlog/3) will not be constant with distance as 
is suggested by the observations. In models for 
which the optical depths are taken to be very 
large, we find that the slope of the R(3) line is the 
largest. As a consequence, the line ratios will vary 
greatly as a function of distance from Mira which 
again is not observed. The line ratios R(l)/R(3) 
and R(2)/R(3) increase the further away from the 
star they are measured. (In a less optically thick 
model, the slopes will not differ greatly between 
the three lines and thus neither will the line ra- 
tios.) In our observations, however, we will not 
be able to trace such a trend, having in mind the 
constraints due to signal-to-noise ratio and the un- 
certainties involved. 

The mass-loss rate determined here, 2.5 x 10^^ M0 yr" 
is in good agreement with the mass-loss rate we 
estimated using an analytic approach to Mira's 
wind (Ryde et al. 2000). Note, however, that in 
this approach the level densities of CO are not 
known and a mean is estimated, which will lead 
to an uncertainty, cf. Fig. 3. Also, an assumption 
in the derivation of the analytical expression for 
the mass-loss rate is a tangentially optically thin 



wind, which is not strictly the case for a spher- 
ically symmetric and homogeneous wind. In a 
tangentially optically thin wind we expect a de- 
cline of the intensity with angular distance, given 
by a power law with an exponent of —3. However, 
even in a tangentially optically thick wind we are 
able to obtain the same dlog I/d\og/3~—3 slope. 

The asymmetries in the radio lines, which 
are optically thick, are caused, in our 'standard 
model', by self absorption along the line-of-sight. 
Thus, we are able to reproduce the intensities and 
radio-line profiles reasonably well without invok- 
ing a double-wind scenario as has previously been 
done in the literature (cf. e.g., Knapp et al., 1998). 

7. CONCLUSIONS 

In order to reproduce the intensity levels and 
the angular dependence of the intensity of the 
vibrational-rotational lines of the CO fundamental 
band at 4.6 /xm (w=1^0), it is necessary to require 
low radial optical-depths to the point of scatter- 
ing. In a standard CSE model, i.e., a smooth and 
spherically symmetric wind with an inner radius 
located fairly close to the central star, the radial 
optical depths are very large. The photons are 
scattered already close to the star and, due to the 
velocity field, shifted out of the absorption profile, 
leading to an escape. The CSE model that is able 
to meet these new observational constraints rea- 
sonably well needs a form of inner cavity, devoid 
of gaseous CO. The stellar light is then able to 
reach further out before it is scattered; the op- 
tical depths are reduced significantly. The two 
CSE models presented here both require a mass- 
loss rate of 2.5x 10^^ M0 yr^^ and a turbulent ve- 
locity of 1.5kms~^ given an expansion velocity of 
the wind of 2 . 5 km s~ ^ . The only parameter that is 
changed between the two models is the location of 
the inner radius. In the standard model, that best 
reproduces the observed radio and FIR emission, 
we use an inner radius located at approximately 
' 13 R, while for the the cavity model, adopted to 
meet the NIR observational constraints, its loca- 
tion is approximately 50 R*. 

We find that the ratio between the turbulent 
velocity and the expansion velocity, is of great im- 
portance. The flux at 2" is sensitive to this ra- 
tio; the lower the ratio the lower the observed flux 
since more photons are lost 'sideways' due to ve- 



locity shifts out of the absorption profile. Further- 
more, we conclude that for modelling vibrational- 
rotational lines the approximation of CRD is not 
a bad assumption as long as the optical depths are 
not too large. 

The NIR observations in combination with a 
model of the CSE is providing us with new pieces 
of information on the inner-most regions of the 
wind. The properties of this inner wind, that is 
simulated in the model by a void of CO, can be 
explained physically by structures of the gas close 
to the star (within approximately 2") that arc in 
such a form that light is able to pass through. Fur- 
ther observations of the gas in the stellar wind of 
Mira slightly closer to the star could therefore be 
rewarding. 

The inner parts of the circumstcUar environ- 
ment are complicated regions of which we have a 
very limited understanding (see e.g., Ryde et al., 
1999a), and where many of the basic assumptions 
of our model are certainly not valid. Thus, the 
model could fail in describing reality. We need fur- 
ther observational constraints. The Phoenix spec- 
trometer makes it possible to use long-slit spec- 
troscopy, which makes spatially resolved spectra 
along the slit attainable. The slit is then used to 
map the CSE. However, no intcrferometric map- 
pings of the CSE of Mira with a high enough spa- 
tial resolution exist. Plancsas et al. (1990b) pre- 
sented observations with a resolution of 7.8" x 5.4" 
using the Owens Valley millimeter-wave interfer- 
ometer. The IRAM interferometer at Plateau de 
Bure can provide a spatial resolution of under 1". 

We are very greatful for inspiration from, and 
enlightening discussions with B. Gustafsson, H. 
Olofsson, and K. Eriksson. The referee is thanked 
for valuable suggestions. This research was sup- 
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Swedish Natural Science Research Council, and 
the Royal Swedish Academy of Sciences. 
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Fig. 1. — The o Ceti light-curve. 1998 October 29 corresponds to a Juhan date of 2451116. At that time 
Ceti had a visual magnitude of 7.8 and was in a pre-maximum phase. (Mattel, J. A., 1999, Observations 
from the AAVSO International Database, private communication). 
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Fig. 2. — Modelling results for the best standard model of Mira's wind, a-c) The JCMT millimeter-wave 
observations, in main-beam brightness temperature, overlayed with the model line-profiles. The beam size 
used is indicated in the upper right-hand corner of each panel. Note that for the C0( J=3^2) we also show 
a scaled line profile for easier comparison with the observed shape of the emission line, d) ISO observations 
and the model overlays using a beam size of 70". 
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Fig. 3. — Modelling results for the best standard model of Mira's wind, a) The intensity of the NIR 
emission lines R(l) (full), R(2) (dashed), and R(3) (dash-dotted) as a function of the distance (in cm and 
arcscconds) from the star, b) The kinetic temperature structure of the model of the circumstellar envelope, 
given by a power law, is shown as a thick full line. The three other lines show the calculated excitation- 
temperature structure of the infrared R(2)-line (dotted), the radio CO(J=3^2) line (dashed), and the 
far-infrared CO(J=18^17) line (dash-dotted). All lines are non-thermal. The CO(J=18^17) line is shown 
only out to r=lxl0^^cm since beyond this radius these levels are scarcely populated, see Fig. 3^. c) The 
tangential optical depths of three selected lines (symbols same as in b) . These give an indication of how far 
out in the wind the lines are formed, d) Relative level-populations according to our best model of Mira's 
wind. In the model we include 30 rotational levels each in the two lowest-lying vibrational states. For clarity, 
only 10 rotational levels from each of the two vibrational levels are shown. The upper set of graphs represent 
the populations of the rotational levels in the ground vibrational state (w=0), i.e., J=0— 10. The lower set 
of graphs represent the excited vibrational state v~l. For instance, the R(l) line is a transition from (w=l, 
J=2) to (u=0, J=l). 
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Table 1 
The parameters adopted for the CSE around o Ceti 



Model parameter 



Value 



T* 


2400 K 




Td 


HOOK 




Lho\ 


8900 Lq 




L, 


7100 Lq 




U 


1800 Lq 




R, 


3.8x1013 cm 




Rd 


3R* 




D 


128 pc 




M 


2.5xlO-^M0yr- 


-1 


/co 


5x10-4 




Vc 


2.5kms-i 




^P 


2.0x1016 cm 




^i, standard 


5.0x1014 cm 




^i, cavity 


2.5x1015 cm 




Vt 


1.5kms-i 




To 


150K 




Ro 


1.0x1015 cm 
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Table 2 

Intensities in 10~^^ergs~^cm~2arcsec~^ (measured at 2") , the line intensity ratios 

(averaged over 2 — 3.4"), and the variation of the R(2) intensity with distance from the 

STAR. 

/[R(l)] /[R(2)] /[R(3)] Igl igl MM^ 

Observed: 2.1±1.0 1.5±0.8 1.2±0.6 1.7±0.6 1.2±0.2 -2.8±0.4 

Modelled ('standard'): 0.24 0.16 0.11 2.2 1.5 -3.4 

Modelled ('cavity'): 1.21 1.18 1.12 1.3 1.1 -3.0 



17 



Table 3 
The integrated radio flux in Kkms^\ and the averaged ISO flux in IO^^^Wcm" 

C0(J=2^1) CO(J=3^2) CO(J=4^3) ISO flux 

Observed: 44±9 

Modelled ('standard'): 50 
Modelled ('cavity'): 72 



94±19 


77±15 


2±1 


67 


71 


2.0 


147 


178 


2.1 
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